Effect of a gap on the decoherence of a qubit 
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We revisit the problem of the decoherence and relaxation of a central spin coupled to a bath 
of conduction electrons. We consider both metallic and semiconducting baths to study the effect 
of a gap in the bath density of states (DOS) on the time evolution of the density matrix of the 
central spin. We use two weak coupling approximation schemes to study the decoherence. At low 
temperatures, though the temperature dependence of the decoherence rate in the case of a metallic 
bath is the same irrespective of the details of the bath, the same is not true for the semiconducting 
bath. We also calculate the relaxation and decoherence rates as a function of external magnetic 
fields applied both on the central spin and the bath. We find that in the presence of the gap, there 
exists a certain regime of fields, for which surprisingly, the metallic bath has lower rates of relaxation 
and decoherence than the semiconducting bath. 

PACS numbers: 03.65Yz, 76.20+q, 71.20Nr 



I. INTRODUCTION 

Solid state spin devices are promising candi- 
dates for quantum computation^ and for quantum 
communication 2 ". The basic building block for quan- 
tum computation is the qubit, which is an effective two 
level system. Various realizations of qubits include spins, 
Josephson junction qubits^ and others involving quantum 
dots^. The efficiency of these systems as qubits/quantum 
computing devices depends crucially on their coupling to 
the numerous environmental degrees of freedom. Typi- 
cally, environments are composed of phonons, electrons, 
nuclear spins or other noise inducing objects and lead 
to dissipation and a consequent loss of coherence of the 
qubit. Understanding the dependence of this decoherence 
and relaxation on the physical parameters which define 
the environment is vital to fabricating good qubits. 

The simplest models used to study the effects of dis- 
sipation induced by the environment is the Caldcira- 
Leggett models, and the spin-boson^ model where a cen- 
tral spin couples to a bath of free bosons. The past 
decade has however, seen a lot of studies on various 
kinds of dissipative environments including interacting 
and non-interacting spin baths and electronic baths£~— . 
In the limit of weak coupling between the central spin i.e., 
qubit and the bath, intra-bath interactions were some- 
times seen to have a mitigating effect on the decoherence 
but this is not generically truei. The effect of fermionic 
environments on single and multiple qubit systems has 
also been studied in&i2r— . The ensuing decoherence also 
depends on the nature of the coupling between the qubit 
and the bath fermions. In general, one expects Marko- 
vian decay of the coherence of the qubit at finite temper- 
atures. However, unlike the case of bosonic baths, where 
the decoherence rate increases indefinitely with temper- 
ature, the rate saturates with temperature for spin and 
fermionic baths. This is however, not the case if there is 



long range order present in the bath or when the bath 
exhibits glassy features&i£. 

In this paper, we revisit the problem of decoherence 
induced by a simple fermionic environment weakly cou- 
pled to a qubit i.e., a central spin. We consider a Kondo 
coupling where the central spin interacts with the local 
spin density of the bath fermions at the central spin site 
and focus on the effect of different metallic and insu- 
lating densities of states, in particular, that of a direct 
gap in the density of states. In the weak coupling limit, 
we use the well known time-convolutionless (TCL) and 
Nakajima-Zwanzig (NZ) techniques^ to study decoher- 
ence. Though the two methods are equivalent for metallic 
density of states (DOS), we see that these two approxima- 
tions are inequivalent for a semi-conducting DOS at low 
enough temperatures raising questions about the validity 
of the weak coupling approximations in certain limits^. 
We find that in general the asymptotic decoherence is 
qualitatively the same for a bath with a metallic DOS, 
though details of the DOS are very relevant for deter- 
mining the intermediate time behavior. Furthermore, we 
also calculate the relaxation and decoherence rates, 71 
and 72 as a function of external magnetic fields applied 
both on the environment and the central spin. Though the 
metallic or semiconducting nature of the bath plays an 
important role at low temperatures, the baths are quali- 
tatively indistinguishable at higher temperatures. 

The paper is organized as follows: in Sec. UH we in- 
troduce the model and the weak coupling formalism used 
to obtain the decoherence and relaxation of the central 
spin. We investigate the case of zero magnetic fields in 
Sec. IIIII Our results for the coherence of a qubit coupled 
to a metallic and a semiconducting bath at zero and finite 
temperatures are presented in Sec. IIII Al and Sec. MI Bl 
In Sec. IIV[ we study the dependence of both the relax- 
ation and decoherence rates on external magnetic fields 
applied on the qubit and bath respectively followed by a 
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FIG. 1. Densities of states p U n{E) (dashed lines), p e i{E) (bold 
lines), paq(E) (dot-dashed lines) as a function of E for a metal 
p > A. 



brief discussion of our results in Sec. |V] 



II. WEAK COUPLING FORMALISM 

We present the model Hamiltonian and the method- 
ology used to study the reduced dynamics of the central 
spin interacting with a bath. The Hamiltonian that de- 
scribes a localized spin-i <x c coupled to a bath of non- 
interacting electrons is given by: 



H —Hs + Hb + Hsb 

H = - hoi + E £kaC ka C kc 



.V 



(1) 



where Hs denotes the intrinsic Hamiltonian of the cen- 
tral spin described by the Pauli matrices <x c subjected 
to a magnetic field h, Hb denotes the bath of non- 
interacting electrons and Hsb describes the Kondo cou- 
pling between the central spin and the bath with V = 



<T a pc p p proportional to the spin density of 



the bath electrons at the origin where the central spin 
is positioned. The coupling constant A is smaller than 
all the other scales of the hamiltonian. Cfc a ,ct are the 
fcrmionic annihilation and creation operators with wave 
number k and Eka = £k + oiH (a = ±1) denotes the 
dispersion of the two spin species of the bath electrons 
subjected to an external magnetic field H. 

To study the evolution (decoherence and relaxation) of 
the central spin, we assume that at time t = 0, the central 
spin is in a pure state \ip) = a|t)+/3|^) and that the bath 
is in thermal equilibrium at temperature T. This leads 
to a factorizable initial density matrix f2 = p s (0) 55 Pb 
where p s (0) = \t/}){ip\ and 
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-H B /T 
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where Z = Tr [cxp (— I3Hb)\ is the bath partition func- 
tion. We use units h = fcg = 1 throughout this paper. 
The time evolution of the reduced density matrix of the 
central spin is given by 



p s (t) = Tr B [e" lHt ne iHt ] 



(3) 



where Tr# is the partial trace over the bath degrees of 
freedom. We rewrite the density matrix ([3]) as a Laplace 



Ps(t) ^ 
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where 77 is a real positive number and £ is the Liouville 
operator corresponding to the total Hamiltonian H, i.e., 
CA = [H, A] for any operator A. The density matrix 
can be decomposed in the basis of Pauli spin operators 
as follows: 



Ps(z) 



E 

a—0,x,y,z 



M a (z)a c c 



(5) 



The decoherence and relaxation are given by the quan- 
tities M a which are in general rather difficult to calcu- 
late. However, for weak central spin- bath coupling i.e., 
A — > 0, we can use perturbation theory in conjunction 
with the projection operator technique to calculate the 
components M<J£. These lead to equations for M a which 
go beyond the standard Markovian master equation. To 
second order in A, the M a (z) satisfy the following matrix 
equation 

zMp(z)-J2h<xM a (z)-J2Zfi a (z)M a {z) = (ff£) (6) 

a a 

where T, a p are the self-energies and only h xy = h* x = 2ih 
are nonzero. In the absence of external magnetic fields 
(h = H = 0), since we have a SU(2) invariant bath (i.e., 
£ka = £k-a), the self energy matrix simplifies with the 
only non-zero entries being the diagonal terms E IX = 
T. yy = Yi zz = S. These diagonal self-energies are related 
to the bath correlation function 



E(*) 



-8i 



die 1 



^e{V x {t)V x ) 



(7) 



where for any operator A, A is defined as A(t) = 
e iH B t Ae -iH B t _ ^ and ^ ^ denotes the thermal ex- 
pectation value, i.e., (A) = Tr(Aps) for any operator A. 
Consequently, M x = M y = M z = M with M being given 
by 



M(z) 



z-E(z) 



(8) 



The exact expressions of h a p and for the general 
case arc detailed in the appendix. 

To obtain M(t) and hence p s (t), we first analyti- 
cally continue S to real frequencies to 

lim E(w + irj) = A(w) - iT(u) (9) 

where A and F are related by the Kramcrs-Kronig rela- 
tion. 
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Note that T is even in w whereas A is odd. For the 
case of the electronic bath studied here, calculating the 
correlation in Eq.©, we find that T is given by 

r(w) = 8ttA 2 [S(w) + S(-u)] (11) 

where S(u>) is the dynamical spin structure factor for the 
bath and given by 

S(u) = / dEp{E)p{E + uj)n{E) [1 - n{E + u)] (12) 



0.08 



In Eq.jnj}, p(E) is the density of states (DOS) for the 
bath and n(E) is the Fermi occupation number. 

The inverse Laplace transform of M(z) directly gives 
the Nakajima-Zwanzig (NZ) approximation for the co- 
herence 

1 f°° 

M NZ (t) = - (iwcos(wi)r(u;) (13) 



where 



f(w) = r(w)/ [(w - A(uj)) 2 + T(uf 



(14) 



On the other hand, if the self energy £(z) is analytic in 
the lower half planed, we obtain a simplified form for the 
coherence which is the well known time convolutionless 
approximation (TCL) 

2 f°° sinM/2) 2 , , , , 
IuMtcl (t) ^ / duo K ' > T(lj) 15 

Note that the NZ approximation corresponds to a master 
equation which is non-local in time i.e., it has a mem- 
ory kernel, whereas TCL approximation corresponds to 
a purely local in time master equation 2 ^. Though the 
TCL is local in time, the time-integration over the mem- 
ory kernel is the characteristic feature which can capture 
the non-Markovian nature of the dynamics. If is 
analytic in the lower half plane, both techniques lead to 
the same asymptotic behavior for the coherence but they 
might yield different results in cases where is not 

analytic. As we shall show later, the two methods con- 
verge for the case of metallic DOS wherein S is analytic 
in the lower half plane, but they predict different results 
for the simple case of a semiconducting DOS, because of 
the presence of weak non-analyticities on the real axis. 

In the presence of external magnetic fields, the self- 
energy matrix E^q, has a far more complicated structure 
than the one described above and requires a full matrix 
inversion to obtain the coherences M a . However, the 
two timescales T\ (rate 71) and T2 (rate 72) which de- 
termine the asymptotic relaxation and decoherence i.e., 
(t \ps{t)\f) ~ exp [-t/Ti] and (f \ps(t)\i) ~ exp [-t/T 2 ] 
are given by the bath correlation functions 
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(16) 
(17) 



with V± = V x ± iVy.ln the rest of the paper, we use 



the dimensionless rates 712 = 7i,2/87rA 
h = H = 0, we have 71 = 72 = 7. 



Note that for 




FIG. 2. (color online) Low frequency behavior of 
r(w)/87rA 2 JV 2 as a function of cu (units of A) for different 
values of the chemical potential: p — (black), p — 0.15 A 
(blue), p — 0.2A (red) and temperature: T = (bold lines), 
T = 0.02A (dot-dashed lines), T = 0.05A (dashed lines) for 
the elliptical density with A = 0.1 A. In the Figure, the tem- 
peratures T — 0, T = 0.02A are indistinguishable for p = OA. 
The inset shows the full r(w) /8tt\ 2 N 2 for a metal at p = 0.2A 
(red) and a semiconductor p = (black) at zero temperature. 



III. ZERO MAGNETIC FIELD 



We now use the weak coupling formalism presented in 
the previous section to calculate the coherence M(t) for 
fermionic baths having both metallic and insulating den- 
sity of states. In order to study the effect of the details of 
the DOS on the asymptotic coherence, we consider three 
densities of states: uniform, square root and elliptical, 
each of which has a gap A and a cut-off A: 



-e(\E\-A)Q(A-\E\) 
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Pun{E) ~2(A~A)- 
p sq [E) = Ny/\E\-AS(\E\ - A)9(A - \E\) 



p el {E) = N' 



(A-Af 



E\ 



(A + A) 



(18) 



where the normalization constants A^ 



and 



l(A-A)- 

N' = t^N(A — A)~ 2 have been chosen so that all densi- 
ties of states have the same normalization. Depending on 
the chemical potential p, one obtains a metallic bath for 
p > A and a semiconducting bath for p < A. Though 
real semiconductors have more complicated DOS, we 
nonetheless expect the simple DOS used in this paper 
to capture the essential features of the asymptotic re- 
laxation and decoherence induced by such baths. The 
gaps in semiconductors are typically of the order of 1 eV: 
1.12 eV for Silicon and 0.66 eV for Germanium 2 ^ and the 
cut-off A is of the order of 10 eV. 
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A. Metallic bath 

Here we assume that the chemical potential /i > A. 
At zero temperature and for small frequencies cj, us- 
ing Eqs. (fTTj) and (|12p . we obtain the generic result 
T(u>) = 2K\u>\ reminiscent of the Ohmic behavior 
seen in the standard spin-boson problem^,. The co- 
efficient K depends on the DOS and we find that 
K un = nX 2 /(A - A) 2 , K sq = 4:n\ 2 N 2 (n - A) and 

K el =47rA 2 iV 2 (|f) 2 (A-A)- 1 ( A i-A)(A-/i). The fact 
that K un is independent of /i is an artifact of the uniform 
DOS. This Ohmic bath leads to the expected power law 
decay of the decoherence M(t) ~ t _4jr /'" where for the 
different DOS, the appropriate value of K should be 
substituted. The asymptotic decoherence induced by 
a metallic bath is qualitatively similar to that induced 
by a bosonic bath at T = 0. This correspondence with 
the spin-boson problem does not hold true at finite 
temperatures or in the presence of magnetic fields as 
we will show later. The full behavior of r(a>) for the 
case of the elliptical DOS is plotted in Fig [2] for various 
temperatures and chemical potentials. 
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FIG. 3. (color online) In Mtlc /8nX N as a function of t 
for a fi = (black), 0.2 (red), 0.21 (turquoise) and T = 
(bold- lines), 0.05 (dashed-lines) for the elliptical DOS with 
A = 0.1. In the inset, we plot lnM/JV 2 for NZ and TCL 
approximations for (T,fi) — (0.1,0.2) and A = 0.05. Time, 
gap, chemical potential and coupling strength are in units of 
A. 



At finite temperature, we see two principal features: 
T(u! = 0) 7^ and T(u) is fully analytic. These features 
are also seen for the square root DOS but one encounters 
non-analyticities for the uniform DOS because of the 
discontinuities in the band structure. However, T is 
expected to be analytic for realistic metallic DOS. For 
a given band filling, the low frequency behavior of T is 
found to be qualitatively the same for both square root 
and elliptical density of states, implying that the asymp- 
totic behavior is indeed similar in both cases. For T ^ 0, 
both approximations TCL and NZ predict a Markovian 



decay of the coherence for times t ^> tx ~ 1/T, with a 
decay rate 7 = T(0) as shown in Fig. [31 The equivalence 
of TCL and NZ at all times is clearly illustrated in the 
inset of Fig. [3] For a metallic DOS, the decoherence 
rate 7 oc T at low T and 7 cx TtanhA/T as T — > 00 
where A is the cut-off. For p un (E) at low temperatures 
the exact result is 7 = 2T(2 - er {tl -^)/ T ). Note that 
the temperature independent proportionality factors are 
dependent on the details of the bath DOS. In Fig. 21 we 
plot the rate 7 for metallic bath [i > A for both square 
root and elliptical DOS. As seen in the case of inter- 
acting spin baths^, 7 increases with temperature and 
saturates to a finite DOS-dependent value proportional 
to J p(E) 2 dE at high temperatures. This saturation 
stems from the fact that the electrons have a spin degree 
of freedom and is different from the case of a bosonic 
bath, where the decoherence rate never saturates. FigJ3] 
also shows that a higher chemical potential results in a 
faster decoherence at all temperatures as there are more 
carriers which can dissipate energy. To summarize, in 
the weak coupling limit, both second order TCL and 
NZ schemes predict that a simple metallic bath results 
in a power law decay of the asymptotic coherence at 
T = and a Markovian decay at T ^ irrespective of 
the nature of the underlying DOS. The corresponding 
exponents and the Markovian rate do not depend very 
much on the qualitative details of the DOS. However, 
we expect the intermediate and other short time regimes 
to be indeed dependent on the details of the DOS. 
To understand whether the presence of a gap in the 




FIG. 4. (color online) Dimensionless rate 7/iV 2 as a func- 
tion of temperature in units of A for a semiconductor fi = 
(dashed lines) and a metal [i — 0.2A (bold lines) for the el- 
liptical (black) and square root (red) DOS with A = 0.1 A. 

DOS below the Fermi energy has any impact whatsoever 
on the decoherence, we plot T in the right panel of Fig. 
M for the elliptic DOS for A = and A / at the 
same value of the filling. At T = 0, though the low 
frequency behavior is the same in both plateau 
like structure is indeed seen in the vicinity of uj ~ 2A 
when A 7^ 0. We observe that for asymptotic times, 
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T=0. |l=0.2A. A=0.1A 
T=0, (1=0.1 A. A=0 
T=0.01A, H=0.2A, A=0.1A 
T=0.01A, H=0.1A, A=0 




FIG. 5. Left panel: In M T Lc/8n A 2 iV 2 as a function of t for 
A = (green) and A = 0.L4 (black) at T = (bold lines) and 
T = 0.0L4 (dashed lines). Right panel: Full r(oj)/87rA 2 iV 2 as 
a function of uj for A = and A = 0.1 A at T = 0. 



M(t) oc Ct~ K where C is some constant which depends 
on the bath parameters and A. We find that though 
the exponent k is the same irrespective of the value of 
A, the coefficient C is indeed dependent on A resulting 
in a faster decoherence in the gapless case cf. top left 
panel of Fig. [SJ At low temperatures, as seen in the 
bottom left panel of Fig. [5j though both DOS have 
an asymptotic Markovian regime with the same decay 
rate, the DOS with a gap has an intermediate power 
law regime for times t& ~ 1/ A < t tT which is not 
present for the gapless metallic bath because t& ~ oo. 
This intermediate regime can be quite large depending 
on the parameters of the problem. To summarize, a 
gap in the DOS below the Fermi energy does serve to 
reduce decoherence as compared to the case of a gapless 
DOS and this feature could probably be exploited in 
experiments. 



B. Semiconductor bath 

We now focus on the semiconductor bath with a gap 
at the fermi level. At zero temperature, T(w) = for 
uj < 2 A and T(u) oc A 2 (w - 2A) 2 for u close to 2A 
for all DOS. The full behavior of T(uj) for the elliptical 
DOS is plotted in Fig. H We first focus on the TCL 
results. At T = 0, we find that due to the presence of 
a gap in T, the central spin only decoheres partially. A 
straightforward evaluation of (fT5|) shows the existence of 
three regimes: a short time regime tA <SC 1/A where the 
central spin initially decoheres as a Gaussian followed 
by an oscillatory regime for t < t& with a well defined 
minimum around t ~ 4/(2A + A) and the asymptotic 
regime for t > t&, where M(t) — > const as a power law. 
These three regimes are expected to be generic and are 
indeed seen in the full numerical results for the coherence 



for both square root and elliptical DOS (see Fig. [3J. 

As temperature increases, the thermal activation of 
the gap results in a T which is no longer gapped (cf. Fig. 
[3]). Within the TCL formalism, this immediately implies 
that the asymptotic behavior is Markovian In Mtcl (i) ~ 
— jt. For low enough temperatures T <C A which is the 
real semiconducting regime, we have the usual Gaussian 
behavior for t < tA, followed by an oscillatory regime for 
<A "C t <C t a , a power-law intermediate regime for t& <C 
t -C tx and finally the asymptotic Markovian regime for 
t > t T (cf. Fig. [J). At high temperatures T > A, there 
is no difference between the metal and the semiconductor 
and the intermediate power law regime ceases to exist. 

In contrast to the metallic case, where the decoher- 
ence rate seemed to have the same qualitative behavior 
for all the DOS considered, the rate in the semiconduct- 
ing case seems to be DOS dependent as T — > 0. In this 
limit, 7 oc T 2 exp — (A — /i)/T for p sq (E) and p e i(E) and 
7 oc Tcxp— (A — fx)/T for p un (E). Comparing with the 
case of the metallic bath where 7 oc T, we see that the 
semiconducting bath has a much smaller rate and hence 
longer coherence times which makes it a better environ- 
ment for qubits than a simple metal. However, as T 
increases, 7 has the same qualitative behavior as indi- 
cated in Fig. [4] At higher temperatures, since the gap is 
completely smeared by thermal effects as expected there 
is no real difference between the metallic and the semi- 
conducting baths. 

We now address the question of whether NZ predicts 
similar results for the asymptotic decoherence in the 
semiconducting case. Substituting A(w) obtained from 
Eq. dTUJ) in ([H]). we find that f(w) typically has a 
three peak structure. At low temperatures, r(w) can be 
rewritten as the sum of two disjoint contributions V(uS) = 
r'(w) + r s (w) where T l is a Lorentzian peak at low fre- 
quencies and T s describes two small satellite peaks which 
exist for uj > 2A. This separation of spectral weight in 
f leads to M NZ (t) = M l NZ (t) + M s NZ {t) where the first 
term is the usual Markovian decay M l NZ (t) ~ exp — "ft 
with 7 being the same as that given by the TCL ap- 
proximation. To estimate the correction Mf jz (t) stem- 
ming from the satellite peaks, we note that for ui ^ 2A, 



A(e) ~ Ao+eAi- 
Consequcntly, 



-£ 2 A 2 +e 2 lne+0(e J ) where s = w— 2A. 



7T 2 A 2 



;(LU-2A)' 



(Ac - 2A) 2 

Using (fH?)) and (JT3J), we find that for long times 
1 



(19) 



M s NZ (t)/7ry 



4A : 



■(E c - 2 A) 2 sin (E c t)/t (20) 



where E c > 2 A is a cut-off scale for which r s (w) has a 
quadratic behavior. The true asymptotic behavior pre- 
dicted by the NZ scheme now depends on the competition 
between the exponential and the power law terms. This 
defines a new time scale t nz oc 1 / 7 \og{\ / B^f) where B is 



the pre- factor of the power law term seen in Eq. (|20[) such 
that for tr < t < t nz one obtains the usual Markovian 
behavior, whereas for t 3> t nz one sees a non- Markovian 
power law decay of the coherence. In the weak coupling 
formalism used here, A is expected to be much smaller 
than all the other scales of the hamiltonian, so a t nz of the 
order of A~ 2 does correspond to very very long times for 
which the coherence has already decayed to infinitesimal 
values and is practically unobservable. A full numerical 
integration of (fT"3|) with the appropriate T does show the 
existence of these two asymptotic regimes Mfj Z {t) and 
M l NZ {t) at low temperatures. For higher temperatures, 
the gap in T is completely smeared by thermal fluctua- 
tions and one recovers the pure Markovian decay. The 
preceding analysis shows that the two approximations 
TCL and NZ do not necessarily predict same results in 
the weak coupling limit even for simple environments like 
the semiconductor considered in this section. This di- 
vergence between the TCL and NZ approximations has 
been seen in other problems like spin baths where exact 
solutions are knowr*^. This shows that one should ex- 
ercise caution when using these approximation schemes 
to calculate coherences in the cases of baths with more 
complicated characteristics. 
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FIG. 6. (color online) Dimensionless relaxation rate 
yi(H)/N 2 as a function of temperature in units of A for 
a semiconductor /j, — 0.05,4 (dashed lines) and a metal 
fi = 0.15,4 (bold lines) for fields H = (black) and H = 0.04A 
(orange) for the elliptical DOS with A = OA A. Inset: Di- 
mensionless decoherence (dashed lines) and relaxation (bold 
lines) rates "fi,2(H)/N 2 as a function of 2H for a semicon- 
ductor /i = 0.05,4 (red) and a metal fj, — 0.15,4 (black) at 
T = 0.05,4. 
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FIG. 7. Dimensionless decoherence and relaxation rates 
71,2(^1) /N 2 for a semiconductor /i = (dashed lines) and a 
metal /x = 0.2,4 (bold lines) as a function of temperature 
for different values of h, H — and an elliptical DOS with 
A = OA A. Femperature is in units of A. 



play between the magnetic fields can be used as a way 
to control decoherence and relaxation. Calculating the 
in-plane correlation functions (V XjV (t)V Xj y) and the per- 
pendicular correlation function (V z {t)V z ) and using them 
in Eqs. (fl"6]l and (fT7|) . the dimensionless decoherence and 
relaxation rates are given by: 
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72 
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= Y. I dE P( E )p( E + 2 P( H + h)) (21) 
P =±i J 

x n(E + pH) [1 - n(E + p(2h + H))] 
\ E / p(E) 2 n(E+pH)[l-n(E + P H)} 

(22) 



p=±i 1 



In the following, we assume h,H > 0. We clearly see 
that for H = h = 0, 71 = 72. For the special case of the 
Ising coupling V = (14,0,0) we have 72 = 71/ 2 for all 
fields. As we shall see below, we have a rich spectrum 
of behavior for the two rates as a function of field for 
all DOS. Though the plots FigJBHH] show the numerical 
results for the elliptical density p e i (E) only, we find that 
the key qualitative features summarized below are true 
for all the DOS considered in this paper. The qualita- 
tive features of these results are seen for all the other 
DOS as well. Additionally, we focus on low and inter- 
mediate temperatures, since the metallic (/i > A) and 
semiconducting (/i < A) baths are indistinguishable at 
high temperatures. 



IV. NON-ZERO MAGNETIC FIELD 

In this section, we study the dependence of the asymp- 
totic decoherence and relaxation rates on magnetic fields 
H and h applied on the bath electrons and the central 
spin respectively cf. Eq[TJ We explore whether an inter- 



A. h = 0,H=£0 

Our results for the two rates 71 and 72 as a function 
of H and T are plotted in Figs. |6][8]for both metallic and 
semiconducting baths. We note that contrary to naive 
expectations, an external magnetic field applied on the 
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bath does not always reduce the rates as compared to the 
H = case and in fact show a lot of anomalous behaviors. 
As T — > 0, both 71,72 — > for all H indicating non- 
Markovian behavior. At finite temperatures, as shown in 
the inset of Fig. |6l a common feature for the metallic and 
semiconducting baths is that both rates 71 (H) and 72 (H) 
have a non-monotonic dependence on H regardless of the 
value of the chemical potential \i. This non-monotonicity 
of the rates stems from the presence of a gap in the DOS 
since at high fields one is sensitive to the detailed nature 
of the DOS cf. Eq. tp?!]) . For example, such features are 
not seen in the case of a banal metallic bath described 
by p(E) oc 6(^4 — \E\). 

At low temperatures, a straightforward calculation 
shows that for the semiconductor, ^yf(H) < 7 s and for 
the metal, 7f (H) < 7" 1 if H < A. Here 7 s and 7™ de- 
note the corresponding H = rates for the semiconduct- 
ing and metallic baths discussed in the previous sections. 
On the other hand though 7f < 7 s and 7™(H) < 7™ 
forO<H < \A-fi\, if \A-[j,\ < H < A then 7f(H) > 7 s 
and 72™(i/) < 7™. (see Table [TJ) . This difference between 
the metallic and semiconducting baths stems from the 
fact that for H > |A — fx\, changing the field is similar 
to changing the chemical potential (see eq. [UJ) and the 
bath switches from a metal-like behavior at low temper- 
atures to a semiconductor-like behavior or viceversa. For 
extremely large fields H > A, the relaxation rate goes to 
zero as anticipated but the dccohcrcncc rate 72 remains 
finite. This is true also for h ^ but the critical field 
for which it happens is H = A — h. Note that in this 
high field limit an Ising like coupling would lead to both 
rates 7^2 = 0. The Markovian decoherence accompanied 
by zero rate of relaxation is a manifestation of the non- 
Ising nature of the coupling between the central spin and 
the bath. To summarize, fields H < A tend to partially 
suppress decoherence and relaxation induced by metal- 
lic baths. However, in the case of semiconducting baths, 
though fields H < A — fi do partially suppress relaxation 
and decoherence, fields | A— /i| < H < A tend to augment 
decoherence. 



B. H = 0,h/0 

In this case, where the central spin has intrinsic dy- 
namics, Eq. (j!?Tj) shows that 71 = T(ui = 2h) and 72 = 
r(2/i)/2 + r(0)/2 (cf. Fig. [2] ). Consequently, both 
7i(/i) and 72(^1) have the same functional dependence 
on the central spin field h. We illustrate some scenar- 
ios where the intrinsic dynamics play an interesting role 
in Fig. [Jj For a metallic bath, turning on a small field 
< 2h < A — fi\ at T = has the dramatic effect of 
transforming the non-Markovian decay of the decoher- 
ence and relaxation seen for h = into a Markovian 
decay since 7™ 2 oc h for small h. At low temperatures 
(see Table H|, for fields < 2h < A - fx, j^ 2 (h) > 7™ 
whereas 7f 2 (^-) < 7 s due to the activation of the gap in 
r cf Fig. [2] Moreover, these rates increase monotonically 



with T and saturate to a finite value as T — > 00. 

The intermediate field regime yields the rather sur- 
prising result that for both metallic and semiconducting 
baths 71 initially decreases with temperature, and then 
increases to attain a constant temperature independent 
value (blue curve Fig. [7]). However, 72 augments con- 
siderably in the same range of temperatures considered 
(red curve Fig. [JJ. Additionally, there exists a special 
value of the intrinsic field h c for which the relaxation 
rate (blue curve Fig [7J has a reduced temperature de- 
pendence. This last feature also survives in the presence 
of a bath field H . For higher values of h, 71 monoton- 
ically reduces with temperature but 72 always increases 
with temperature (yellow and purple curves Fig[7J). For 
all values of h, we have 71 > 72 for T < A and 72 > 71 
for higher temperatures. To our knowledge, the possibil- 
ity of an almost temperature independent relaxation at 
intermediate fields and the scenario where the relaxation 
rate decreases and the decoherence rate increases with 
temperature have not been discussed in the literature. 



Interval/bath Metallic Semiconducting 

< H < |p - A\,h = 7^ 2 (H) < 7^ jl^jH) < % 2 
\ti - A| < H < A,h = 7 ™(ff) < 7™ 7| (H) < 7 f, 7 |(H) > % 
< 2h < \fi ~ A I , H = 1 7j" 2 (h) > 7^3 I 7 f, 2 (h) < ^ 

TABLE I. Comparison between decoherence and relaxation 
rates 'y™^ (H, h) and y'l'-f f° r l° w temperatures and only one 
applied field. 
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FIG. 8. (color online) Dimensionless relaxation (bold lines) 
and decoherence rate (dashed lines) ^i^/N 2 as functions of 
2h for a semiconductor bath /i = 0.05A (left panel) and a 
metallic bath /x = 0.15^4 (right panel) for T = (black) and 
T = 0.1A (blue) for the elliptical DOS with A = 0.1A and 
H = (top panel) and H = 0.08A (bottom panel). 
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FIG. 9. (color online) Dimensionless relaxation (bold lines) 
and decoherence rate (dashed lines) 71,2 /N 2 as functions of 
T for a semiconductor bath fj, — 0.05^4 (left panel) and a 
metallic bath \i = 0.15j4 (right panel) for H = (black) and 
H = 0.08,4 (red) for the elliptical DOS with A = 0.1A and 
2h = 0.1A. 



C. h^0,H^0 

In the presence of both fields, the low temperature 
behavior of the rates is an amalgam of the two cases 
studied above and is sensitive to the nature of the 
bath. The results for the relaxation rate ji(H,h) are 
plotted in Fig. [5] and Fig. |H] for both baths. For a 
given value of H and T, the decoherence and relaxation 
rates have the same functional dependence on h. We 
now consider the T = case. Firstly, we note that 
at T = 0, ji(H,h) > 72 (if, ft.) for all reasonable fields 
irrespective of the nature of the bath. Moreover, for 
the metallic bath, if < H < |A — fi\ then both 7^2 
are non-zero for any finite field h. If the field is further 
increased |A — < H < A + /1, then 7^2 = for 
2h < fi + A — H. This is clearly illustrated in the right 
panels of Fig. [8] As a result, the relaxation will be 
Markovian or non-Markovian depending on the values 
of h and H. The case of the semiconducting bath 
is the opposite: 71.2 show a gap as a function of the 
field h for all H < A + fx, indicating non-Markovian 
behavior at low h and Markovian behavior for fields 
2/i > A + — H as shown in the left panels of Fig. [3] 
Finally, for large fields H > fi + A any finite h leads to 
finite rates irrespective of the value of [i. Consequently, 
we have the surprising result that there exists a range 
of fields max(A + fi s — H, 0) < 2h < A + /i TO — H and 
A — fi m I < H < A + ji m where fi s and fi m denote the 
corresponding chemical potentials for the semiconduct- 
ing and metallic baths, for which the metallic bath is 
more effective in suppressing relaxation and decoherence 
than the semiconducting bath. This feature persists at 
low temperatures T -C A + fx s — H, wherein there is a 
regime h < h±(T), where the semiconductor bath has 
lower rates as expected and a regime hi(T) < h < h-2.iT), 
where the metallic bath has lower rates as compared 
to the semiconductor. The interval — hi — > as 
temperature increases. In the first regime we find that 
72 > 71 for all baths. We also see from Fig. [5] that 
for < H,h < |A - n\ ^f(H,h) < 7l s (iJ) < 7 s and 



ji(H,h) < 7 s . For larger fields |A - fi\ < H < A 
and small enough h, 7f(-H, h) > 7 s reminiscent of 
h = case. For a metallic bath, if < H < A — fi\ 
and 2h < \i — A — H then both the decoherence and 
the relaxation rate are greater than the decoherence 
and relaxation when no fields are present and we 
have 7^ 2 (h) > ( H > h ) > 7™- At slightly higher 
temperatures, the semiconducting bath has a lower 
rate and at very high temperatures the two baths are 
indistinguishable as expected. 

To summarize, we find that external magnetic fields 
do not necessarily reduce the relaxation and decoherence 
rates of the qubit. Depending on the nature of the bath 
and the values of the field, one sees a rich and varied 
behavior of the rates. For instance, in the absence of an 
intrinsic qubit field (h = 0), we find that the presence of 
a gap in the density of states leads to a non-monotonic 
variation of the rates as a function of the bath field H . 
In the presence of both fields, h and H, an interplay be- 
tween the two magnetic fields leads to a very interesting 
regime wherein for a qubit field h\ < h < hi, the metallic 
bath has lower rates as compared to the semiconductor. 
The bounds h\ and hi are functions of the gap A, the 
bath field H and the temperature. For any given H , 
there exists a value of the bath field h c for which we have 
an almost temperature independent relaxation. For very 
high bath fields, the relaxation rate goes to zero but the 
decoherence rate remains finite signaling the non-Ising 
nature of the qubit-bath coupling. It would be very in- 
teresting to see if these features survive in the presence 
of an additional transverse field on the qubit. 



V. DISCUSSION 

We have studied the decoherence and relaxation of a 
central spin weakly coupled to a bath of electrons de- 
scribed by different densities of states having a gap, using 
two well known approximation schemes, NZ and TCL. 
Though both methods predict an asymptotic Markovian 
decay for a metallic bath, we have shown that in the case 
of a semiconducting bath, the TCL predicts a Markovian 
behavior at low T whereas the NZ approximation pre- 
dicts non-Markovian behavior. The study of the validity 
of both approximation techniques is left for future work. 
We emphasize that one has to exert care when using ap- 
proximation schemes to study the time evolution of the 
density matrix. We find that a gap in the spectrum, even 
one away from the Fermi level, generates an intermedi- 
ate power law regime for the coherences at finite tem- 
peratures. Depending on the parameters of the problem, 
this power law regime can be made sufficiently large so 
as to be relevant for experiments. One might also ques- 
tion the validity of the weak coupling formalism used here 
since it is well known that even the equilibrium physics of 
the Kondo model is governed by strong coupling physics, 
which leads to a complete screening of the Kondo spin. 
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This was explored in Ref£, where numerical renormaliza- 
tion group methods were used to study the decoherence 
of the Kondo spin. It was found that for any value of the 
coupling A there is a Kondo timescale t\ oc exp const /A 
beyond which the Kondo spin gets screened dynamically. 
The weak coupling results discussed here were indeed 
valid for intermediate times t <C t\. In this paper, since 
we are only interested in very small A, we expect our re- 
sults to be valid for realistically long times since the time 
scale t\ —> oo as A — > 0. 

We have also calculated the relaxation and decoherence 
rate as a function of external magnetic fields applied on 
the bath and/or central spin. We encounter novel sit- 
uations where for moderately large fields on the central 
spin alone result in a nearly temperature independent re- 
laxation but a rapidly increasing decoherence. We find 
that the presence of gap anywhere in the density of states 
has important ramifications for the rates. This leads to 
the case where for a certain range of the bath field H, 
there is a regime h\(T) < h < fi2(T), where we have the 
surprising result that the metallic bath has lower rates of 
relaxation and decoherence as compared to the semicon- 



The hp a are given b 

h l3a = ^ 

Where Tr it the total trace. For the specific case of the 
Kondo coupling: 

hp a = -Tr(crp[Hs,(T a ]) (A3) 
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Appendix A: Self-energy 
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